Chapter 9 

Atomic structure 



Previously, we have seen that the quantum mechanics of atomic hydrogen, and 
hydrogen-hke atoms is characterized by a large degeneracy with eigenvalues 
separating into multiplets of n^-fold degeneracy, where n denotes the principle 
quantum number. However, although the idealized Schrodinger Hamiltonian, 

Ho^^ + V{r), V{r) = --^—, (9.1) 
2m 47reo r 

provides a useful platform from which develop our intuition, there are several 
important effects which mean that the formulation is a little too naive. These 
"corrections" , which derive from several sources, are important as they lead to 
physical ramifications which extend beyond the realm of atomic physics. Here 
we outline some of the effects which need to be taken into account even for 
atomic hydrogen, before moving on to discuss the quantum physics of multi- 
electron atoms. In broad terms, the effects to be considered can be grouped 
into those caused by the internal properties of the nucleus, and those which 
derive from relativistic corrections. 

To orient our discussion, it will be helpful to summarize some of key aspects 
of the solutions of the non-relativisitic Schrodinger equation, Hqi/; — Etp on 
which we will draw: 

Hydrogen atom revisited: 

> As with any centrally symmetric potential, the solutions of the Schrodinger 
equation take the form V'^m(r) = R{r)Yg^{6, (j)), where the spherical har- 
monic functions Y^miO, (f>) depend only on spherical polar coordinates, 
and R{r) represents the radial component of the wavefunction. Solving 
the radial wave equation introduces a radial quantum number, > 0. 
In the case of a Coulomb potential, the energy depends on the principal 
quantum number n = rij. + i + 1 > 1, and not on and i separately. 

The fine structure constant 
known to great accuracy and 
given by, 

a = 7.297352570(5) x 10" 
1 

~ 137.035999070(9) ' 



> For atomic hydrogen [Z — 1), the energy levels of the Hamiltonian (9.1) 
are given by 

^ Ry ( \ m e^l I22 

where ao — ^p^^ is the Bohr radius, a = jf^^ denotes the fine 
structure constant, and strictly speaking m represents the reduced 
mass of the electron and proton. Applied to single electron ions with 

higher atomic weight, such as He"*", Li^"*", etc., the Bohr radius is reduced 
by a factor 1/Z, where Z denotes the nuclear charge, and the energy is 
given by = -jRy = -^mc^{Za)'^. 
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> Since n > 1 and > 0, the allowed combinations of quantum num- 
bers are shown on the right, where we have introduced the conventional 
notation whereby values of ^ = 0, 1, 2, 3, 4 • • • are represented by letters 
s,p,d, f, g ■ ■ ■ respectively. 

£> Since En depends only on n, this implies, for example, an exact degen- 
eracy of the 2s and 2p, and of the 3s, 3p and 3d levels. 

These results emerge from a treatment of the hydrogen atom which is 
inherently non-relativistic. In fact, as we will see later in our discussion of the 
Dirac equation in chapter 15, the Hamiltonian (9.1) represents only the leading 
term in an expansion in v^jc? ^ (Za)^ of the full relativistic Hamiltonian 
(see below). Higher order terms provide relativistic corrections, which impact 
significantly in atomic and condensed matter physics, and lead to a lifting of 
the degeneracy. In the following we will discuss and obtain the heirarchy of 
leading relativistic corrections.^ This discussion will provide a platform to 
describe multi-electron atoms. 



9.1 The "real" hydrogen atom 

The relativistic corrections (sometimes known as the fine-structure cor- 
rections) to the spectrum of hydrogen-like atoms derive from three different 
sources: 
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To see that u^/c^ {Zaf , we 
may invoke the virial theorem. 
The latter shows that the aver- 
age kinetic energy is related to 
the potential energy as (T) = 
— Therefore, the aver- 

age energy is given by {E) = 
(T) + (V) = -(T) = -im«2. 
We therefore have that = 
Ry = ^mc^{Za)^ from which 
follows the relation v^/c^ ~ 
(Za)2. 



i> relativistic corrections to the kinetic energy; 

i> coupling between spin and orbital degrees of freedom; 

> and a contribution known as the Darwin term. 

In the following, we will discuss each of these corrections in turn. 

9.1.1 Relativistic correction to the kinetic energy 

Previously, we have taken the kinetic energy to have the familiar non-relativistic 
form, However, from the expression for the relativistic energy- momentum 
invariant, we can already anticipate that the leading correction to the non- 
relativistic Hamiltonian appears at order p^. 



2^2 



r^-^-^ — 2^ 2 , P 1 (P ) 

V p + m^c^ — mc H 5-„ 

2m 8 m'^c^ 



As a result, we can infer the following perturbation to the kinetic energy of 
the electron, 



I(P^) 



2\2 



8 m^c^ 



When compared with the non-relativistic kinetic energy, p^/2r 
that the perturbation is smaller by a factor of p^ /rn?c^ = jc 
Hi is only a small perturbation for small atomic number, Z <C 1/a ~ 137. 



one can see 



I.e. 



^It may seem odd to discuss relativistic corrections in advance of the Dirac equation 
and the relativistic formulation of quantum mechanics. However, such a discussion would 
present a lengthy and unnecessarily complex digression which would not lead to further 
illumination. We will therefore follow the normal practice of discussing relativistic corrections 
as perturbations to the familiar non-relativistic theory. 
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We can therefore make use of a perturbative analysis to estimate the scale of 
the correction. 

In principle, the large-scale degeneracy of the hydrogen atom would de- 
mand an analysis based on the degenerate perturbation theory. However, 
fortunately, since the off-diagonal matrix elements vanish,^ 

{nem\Hi\ni'm') =0 for £ ^ i' or m m' , 

degenerate states are uncoupled and such an approach is unnecessary. Then 



making use of the identity, Hi 



2mc' 



[Ho — V{r)] , the scale of the resulting 



energy shift can be obtained from first order perturbation theory, 

{HiUm ^ {nim\Hi\nim) - {E^ - 2E^{V{r))ni + {V\r))nif • 

Since the calculation of the resulting expectation values is not particularly 
illuminating, we refer to the literature for a detailed exposition^ and present 
here only the required identities (right). From these considerations, we obtain 
the following expression for the first order energy shift. 



{Hi)nlm — 



mc 



Za 
n 



n 

'+1/2 



(9.2) 



By making use of the form of the 
radial wavefunction for the hy- 
drogen atom, one may obtain the 
identities, 



aon^ 



nt 



'+1/2) 



From this term alone, we expect the degerenacy between states with different 
values of total angular momentum I to be lifted. However, as well will see, 
this conclusion is a little hasty. We need to gather all terms of the same 
order of perturbation theory before we can reach a definite conclusion. We 
can, however, confirm that (as expected) the scale of the correction is of order 

^^g"*^"* ~ ('^)^- '^o^ tmn to the second important class of corrections. 
9.1.2 Spin-orbit coupling 

As well as revealing the existence of an internal spin degree of freedom, the 
relativistic formulation of quantum mechanics shows that there is a further 
relativistic correction to the Schrodinger operator which involves a coupling 
between the spin and orbital degrees of freedom. For a general potential V{r), 
this spin-orbit coupling takes the form. 



H2^^\drV^r))S-t. 



For a hydrogen-like atom, V{r) — 



1 Ze' 
47ren r 



, and 



H2 



1 1 



2'm?c^ 47reo 



S L. 



O Info. Physically, the origin of the spin-orbit interaction can be understoon 

from the following considerations. As the electron is moving through the electric 
field of the nucleus then, in its rest frame, it will experience this as a magnetic 
field. There will be an additional energy term in the Hamiltonian associated with 
the orientation of the spin magnetic moment with respect to this field. We can make 
an estimate of the spin-orbit interaction energy as follows: If we have a central field 
determined by an electrostatic potential V{r), the corresponding electric field is given 

^The proof runs as follows: Since [//i, L^] = 0, f {(!((! + 1) - l{l + 1)] {nlm\H^\n(! rW) = 
0. Similarly, since [H-i_,L^\ = 0, fi(m' — m){nim\Hi\n(.' m!) = 0. 
^see, e.g., Ref [1]. 
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by E = — VV(r) = —er{drV). For an electron moving at velocity v, this translates 

to an effective magnetic field B = x E. The magnetic moment of the electron 
associated with its spin is equal to /itg = 5s 2^ S = — ^S, and thus the interaction 
energy is given by 

-M, • B = • (V X E) = -^^S • (p X BridrV)) = ^^-{drV) L • S , 

rnc-' [rricY [mcY r 

where we have used the relation pxer = px^ = — ^. In fact this isn't quite correct; 
there is a relativistic effect connected with the precession of axes under rotation, called 
Thomas precession which multiplies the formula by a further factor of \ . 



Once again, we can estimate the effect of spin-orbit coupling by treating H2 
as a perturbation. In the absence of spin-orbit interaction, one may express 
the eigenstates of hydrogen-like atoms in the basis states of the mutually 
commuting operators, Hq, L^, L^, S^, and Sz- However, in the presence 
of spin-orbit coupling, the total Hamiltonian no longer commutes with 
or Sz (exercise). It is therefore helpful to make use of the degeneracy of 
the unperturbed Hamiltonian to switch to a new basis in which the angular 
momentum components of the perturbed system are diagonal. This can be 
achieved by turning to the basis of eigenstates of the operators, Hq, J^, J^, 
L^, and S^, where J = L + S denotes the total angular momentum. (For a 
discussion of the form of these basis states, we refer back to section 6.4.2.) 

Making use of the relation, = L^ + S^ + 2L-S, in this basis, it follows 
that, 

L-S = ^(P-L^-S^). 

Combining the spin and angular momentum, the total angular momentum 
takes values j — l± 1/2. The corresponding basis states \j = i di 1/2, mj,i) 
(with s — 1/2 implicit) therefore diagonalize the operator, 



S-L\j = £±l/2,mj,e) 



-e-1 



\£±l/2,mj,i) , 



where the brackets index j = £ + 1/2 (top) and j = i — 1/2 (bottom). As for 
the radial dependence of the perturbation, once again, the off-diagonal matrix 
elements vanish circumventing the need to invoke degenerate perturbation 
theory. As a result, at first order in perturbation theory, one obtains 

, „ > 1 / i \ Ze"^ / 1 



1 2 ^-2 



-e-l J 47reo 
Then making use of the identity (right), ^ one obtains 



rd 



1 



2/n,j=e±l/2,mj,e 
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1 



Note that, for £ = 0, there is no orbital angular momentum with which to 
couple! Then, if we rewrite the expression for {Hi) (9.2) in the new basis, 

4 / 1 \ 



n,i=t±\/2,i 



1 



— mc 



Za 
n 



± 



and combining both of these expressions, for £ > 0, we obtain 

{Hi + H2)n,j=l±l/2 



1 2fZaV 

2 \n ) 



while for £ ■ 



j + 1/2 

0, we retain just the kinetic energy term (9.2). 



For details see, e.g., Ref [1]. 



For £ > 0, 



nt 



f mcaZ\ 1 
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9.1.3 Darwin term 

The final contribution to the Hamiltonian from relativistic effects is known as 
the Darwin term and arises from the "Zitterbewegung" of the electron - 
trembling motion - which smears the effective potential felt by the electron. 
Such effects lead to a perturbation of the form, 



/fa 



-Q 



nuclear 



(r) 



47reo 8(mc) 



where Qnuciear(r) = Ze5^^\r) denotes the nuclear charge density. Since the 
perturbation acts only at the origin, it effects only states with £ = 0. As a 
result, one finds that 



47reo 8(mc)' 



47r|V^„(0)|^ 



-mc 



n 



Intriguingly, this term is formally identical to that which would be obtained 
from {H2) dX t — Q. As a result, combining all three contributions, the total 
energy shift is given simply by 



n 

j + 1/2. 



(9.3) 



a result that is independent of (. and nij. 

To discuss the predicted energy shifts for particular states, it is helpful to 
introduce some nomenclature from atomic physics. For a state with principle 
quantum number n, total spin s, orbital angular momentum I, and total an- 
gular momentum j, one may use spectroscopic notation v?^^^Lj to define 
the state. For a hydrogen-like atom, with just a single electron, 2s + 1 = 2. 
In this case, the factor 2s + 1 is often just dropped for brevity. 

If we apply our perturbative expression for the relativistic corrections (9.3), 
how do we expect the levels to shift for hydrogen-like atoms? As we have seen, 
for the non-relativistic Hamiltonian, each state of given n exhibits a 2n^-fold 
degeneracy. For a given multiplet specified by n, the relativistic corrections 
depend only on j and n. For n = 1, we have I — Q and j = 1/2: Both IS1/2 
states, with rrij — 1/2 and —1/2, experience a negative energy shift by an 



amount AEi 



l/2,mj,0 



■\Z'^c? Ry. For n = 2, ^ can take the values of or 1. 



With 2 — 1/2, both the former 25*1/2 state, and the latter 2Pi/2 states share 

Ry, 



2,l/2,mj,l 
1 74„,2 1 



64' 



the same negative shift in energy, A.B2,i/2,mj ,o 

while the 2P3/2 experiences a shift of Ai?2,3/2,mj,i = ~'^^'^oi^^y- Finally, for 
n — i, I can take values of 0, 1 or 2. Here, the pairs of states 351/2 3Pi/2, 
and 3P3/2 and 2D-^j2 each remain degenerate while the state 3D5/2 is unique. 
These predicted shifts are summarized in Figure 9.1. 

This completes our discussion of the relativistic corrections which develop 
from the treatment of the Dirac theory for the hydrogen atom. However, 
this does not complete our discription of the "real" hydrogen atom. Indeed, 
there are further corrections which derive from quantum electrodynamics and 
nuclear effects which we now turn to address. 

9.1.4 Lamb shift 

According to the perturbation theory above, the relativistic corrections which 
follow from the Dirac theory for hydrogen leave the 23^ 12 aiid 2Pi/2 states 
degenerate. However, in 1947, a careful experimental study by Willis Lamb 



Willis Eugene Lamb, 1913-2008 

A physicist who 
won the Nobel 
Prize in Physics 
in 1955 "for his 
discoveries con- 
cerning the fine 
structure of the 
hydrogen spectrum". Lamb and 
Polykarp Kusch were able to pre- 
cisely determine certain electromag- 
netic properties of the electron. 
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3P3,,3D3„ 3P3,„3D3,; Figure 9.1: Figure sliowing tlie 

3D5/2 3D5,; iieirarcliy of energy sliifts of tlie spec- 

3Si,7. 3P„2 3Pi,; I — — tra of liydrogen-liice atoms as a result 

2P3,; 1 — — of relativistic corrections. Tiie first 

2Si,2, 2P,yj ^^^^J column shows the energy spectrum 

2S1/2 f = i predicted by the (non-relativistic) 

Bohr theory. The second column 
shows the predicted energy shifts 
from relativistic corrections arising 
from the Dirac theory. The third col- 
umn includes corrections due quan- 
tum electrodynamics and the fourth 
column includes terms for coupling 
to the nuclear spin degrees of free- 
dom. The H-a line, particularly 
important in the astronomy, corre- 
sponds to the transition between the 
levels with n = 2 and n = 3. 



ISi 



ISin 



Bohr Dirac Lamb fiyperfine 

levels fine structure shift structure 

(increased) (increased) (increased) 



and Robert Retherford discovered that this was not in fact the case:^ 
state is slightly lower in energy than the 2Si/2 state resulting in a small shift 
of the corresponding spectral line - the Lamb shift. It might seem that such a 
tiny effect would be deemed insignificant, but in this case, the observed shift 
(which was explained by Hans Bethe in the same year) provided considerable 
insight into quantum electrodynamics. 

In quantum electrodynamics, a quantized radiation field has a zero-point 
energy equivalent to the mean-square electric field so that even in a vacuum 
there are fiuctuations. These fiuctuations cause an electron to execute an 
oscillatory motion and its charge is therefore smeared. If the electron is bound 
by a non- uniform electric field (as in hydrogen), it experiences a different 
potential from that appropriate to its mean position. Hence the atomic levels 
are shifted. In hydrogen-like atoms, the smearing occurs over a length scale. 



2a 

TT 



h 

mc 
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some five orders of magnitude smaller than the Bohr radius. This causes the 
electron spin (/-factor to be slightly different from 2, 




Hans Albrecht Bethe 1906-2005 

A German- 
American 
physicist, and 
Nobel laureate 
in physics "for 
his work on the 
theory of stellar 
nucleosynthesis." 
A versatile theo- 
retical physicist, 
Bethe also made important contribu- 
tions to quantum electrodynamics, 
nuclear physics, solid-state physics 
and particle astrophysics. During 
World War II, he was head of the 
Theoretical Division at the secret 
Los Alamos laboratory developing 
the first atomic bombs. There he 
played a key role in calculating the 
critical mass of the weapons, and 
did theoretical work on the implosion 
method used in both the Trinity test 
and the "Fat Man" weapon dropped 
on Nagasaki. 
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There is also a slight weakening of the force on the electron when it is very 
close to the nucleus, causing the 2Si/2 electron (which has penetrated all the 
way to the nucleus) to be slightly higher in energy than the 2Pi/2 electron. 
Taking into account these corrections, one obtains a positive energy shift 



amb 



na Ry x — a In — — 
Stt aZ 



for states with i = 0. 

^W. E. Lamb and R. C. Retherfod, 
Method, Phys. Rev. 72, 241 (1947). 



Fine Structure of the Hydrogen Atom by a Microwave 




5 10 GHz iv 

Hydrogen fine structure and fiy- 
perfine structure for the n = 3 to 
n = 2 transition (see Fig. 9.1). 
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9.1.5 Hyperflne structure 

So far, we have considered the nucleus as simply a massive point charge respon- 
sible for the large electrostatic interaction with the charged electrons which 
surround it. However, the nucleus has a spin angular momentum which is 
associated with a further set of hyperfine corrections to the atomic spec- 
tra of atoms. As with electrons, the protons and neutrons that make up a 
nucleus are fermions, each with intrinsic spin 1/2. This means that a nucleus 
will have some total nuclear spin which is labelled by the quantum number, 
/. The latter leads to a nuclear magnetic moment, 

Ze 

^^^^^2il^^- 

where Mjv denotes the mass of the nucleus, and gjv denotes the gyromagnetic 
ratio. Since the nucleus has internal structure, the nuclear gyromagnetic ratio 
is not simply 2 as it (nearly) is for the electron. For the proton, the sole 
nuclear constituent of atomic hydrogen, gp k, 5.56. Even though the neutron 
is charge neutral, its gyromagnetic ratio is about —3.83. (The consitituent 
quarks have gyromagnetic ratios of 2 (plus corrections) like the electron but 
the problem is complicated by the strong interactions which make it hard to 
define a quark's mass.) We can compute (to some accuracy) the gyromagnetic 
ratio of nuclei from that of protons and neutrons as we can compute the 
proton's gyromagnetic ratio from its quark constituents. Since the nuclear 
mass is several orders of magnitude higher than that of the electron, the nuclear 
magnetic moment provides only a small perturbation. 

According to classical electromagnetism, the magnetic moment generates 
a magnetic field 



47rr3 3 

To explore the effect of this field, let us consider just the s-electrons, i.e. ^ = 0, 
for simplicity.^ In this case, the interaction of the magnetic moment of the 
electrons with the field generated by the nucleus, gives rise to the hyperfine 
interaction, 

mc 

For the l — ^ state, the first contribution to B vanishes while second leads to 
the first order correction, 

2^ 8 ml 

N ■ 



{H^y,)n,i/2,o = - na Ry X -^9nj^j-,S ■ I 



Once again, to evaluate the expectation values on the spin degrees of free- 
dom, it is convenient to define the total spin F = I + S. We then have 

pS • I = ^(f2 - S2 - I2) = i(F(F + 1) - 3/4 - /(/ + 1)) 



_ 1 r / F=i 

~2\-/-l F^I 



1/2 
1/2 



Therefore, the Is state of Hydrogen is split into two, corresponding to the 
two possible values F = and 1. The transition between these two levels has 
frequency 1420 Hz, or wavelength 21 cm, so lies in the radio waveband. It 



^For a full discussion of the influence of the orbital angular momentum, we refer to Ref. 
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is an important transition for radio astronomy. A further contribution to the 
hyperfine structure arises if the nuclear shape is not spherical thus distorting 
the Coulomb potential; this occurs for deuterium and for many other nuclei. 

Finally, before leaving this section, we should note that the nucleus is not 
point-like but has a small size. The effect of finite nuclear size can be 
estimated perturbatively. In doing so, one finds that the s {£ = 0) levels are 
those most effected, because these have the largest probability of finding the 
electron close to the nucleus; but the effect is still very small in hydrogen. It 
can be significant, however, in atoms of high nuclear charge Z, or for muonic 
atoms. 



A muon is a particle somewhat 
like an electron, but about 200 
times heavier. If a muon is 
captured by an atom, the cor- 
responding Bohr radius is 200 
times smaller, thus enhancing 
the nuclear size effect. 



This completes our discussion of the "one-electron" theory. We now turn 
to consider the properties of multi-electron atoms. 



9.2 Multi-electron atoms 

To address the electronic structure of a multi-electron atom, we might begin 
with the hydrogenic energy levels for an atom of nuclear charge Z, and start 
filling the lowest levels with electrons, accounting for the exclusion principle. 
The degeneracy for quantum numbers {n,i) is 2 x {2i + 1), where {2£ + 1) is 

the number of available mg values, and the factor of 2 accounts for the spin 
degeneracy. Hence, the number of electrons accommodated in shell, n, would 
be 2 X n^, 

n i Degeneracy in shell Cumulative total 

^1 2 2 

2 0,1 (l + 3)x2 = 8 10 

3 0,1,2 (1 + 3 + 5) X 2 = 18 28 

4 0,1,2,3 (1 + 3 + 5 + 7) X 2 = 32 60 

We would therefore expect that atoms containing 2, 10, 28 or 60 electrons 
would be especially stable, and that in atoms containing one more electron 
than this, the outermost electron would be less tightly bound. In fact, if we 
look at data (Fig. 9.2) recording the first ionization energy of atoms, i.e. the 
minimum energy needed to remove one electron, we find that the noble gases, 
having Z — 2, 10, 18, 36 • • • are especially tightly bound, and the elements 
containing one more electron, the alkali metals, are significantly less tightly 
bound. 

The reason for the failure of this simple-minded approach is fairly obvious - 
we have neglected the repulsion between electrons. In fact, the first ionization 
energies of atoms show a relatively weak dependence on Z: this tells us that the 
outermost electrons are almost completely shielded from the nuclear charge.^ 
Indeed, when we treated the Helium atom as an example of the variational 
method in chapter 7, we found that the effect of electron-electron repulsion 
was sizeable, and really too large to be treated accurately by perturbation 
theory. 

'^In fact, the shielding is not completely perfect. For a given energy shell, the effective 
nuclear charge varies for an atomic number Z as Z^a ~ (1 + a) where a > characterizes 
the ineffectiveness of screening. This implies that the ionization energy Iz = —Ez ~ Z^ff ~ 
(1 + 2aZ). The near-linear dependence of Iz on Z is reflected in Fig. 9.2. 
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Figure 9.2: Ionization energies of tlie elements. 
9.2.1 Central field approximation 

Leaving aside for now the influence of spin or relativistic effects, the Hamilto- 
nian for a multi-electron atom can be written as 



H =y V- - + > 



2m 47reo 



1 e 



2 



. 47reo Tij 

i<j 



where rij = \Yi — Yj\. The first term represents the "single-particle" contri- 
bution to the Hamiltonian arising from interaction of each electron with the 
nucleus, while the last term represents the mutual Coulomb interaction be- 
tween the constituent electrons. It is this latter term that makes the generic 
problem "many-body" in character and therefore very complicated. Yet, as we 
have already seen in the perturbative analysis of the excited states of atomic 
Helium, this term can have important physical consequences both on the over- 
all energy of the problem and on the associated spin structure of the states. 

The central field approximation is based upon the observation that 
the electron interaction term contains a large central (spherically symmetric) 
component arising from the "core electrons" . From the following relation, 

J2 0)1' = const. 

m=-l 

it is apparent that a closed shell has an electron density distribution which 
is isotropic (independent of 6 and cf)) . We can therefore develop a perturbative 
scheme by setting H = Hq + Hi, where 



■^V2- + r, 

Am 47reo 



Here the one-electron potential, Ui{r), which is assumed central (see below), 
incorporates the "average" effect of the other electrons. Before discussing how 
to choose the potentials Ui{r), let us note that Hq is separable into a sum of 
terms for each electron, so that the total wavefunction can be factorized into 
components for each electron. The basic idea is first to solve the Schrodinger 
equation using Hq, and then to treat Hi as a small perturbation. 

On general grounds, since the Hamiltonian Hq continues to commute with 
the angular momentum operator, [iiToiL] = 0, we can see that the eigenfunc- 
tions of Hq will be characterized by quantum numbers (n, £, mi, m^). However, 
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since the effective potential is no longer Coulomb-like, the I values for a given 
n need not be degenerate. Of course, the difhcult part of this procedure is 
to estimate Ui{r); the potential energy experienced by each electron depends 
on the wave function of all the other electrons, which is only known after the 
Schrodinger equation has been solved. This suggests that an iterative approach 
to solving the problem will be required. 

To understand how the potentials Ui{r) can be estimated - the self- 
consistent field method - it is instructive to consider a variational ap- 
proach due originally to Hartree. If electrons are considered independent, the 
wavefunction can be factorized into the product state, 

*({ri}) = Vn(ri)V'i2(i'2) ■ --ipij^irN) , 

where the quantum numbers, ik = {nimgms)k, indicate the individual state 
occupancies. Note that this product state is not a properly antisymmetrized 
Slater determinant - the exclusion principle is taken into account only in so 
far as the energy of the ground state is taken to be the lowest that is consistent 
with the assignment of different quantum numbers, nim^ms to each electron. 
Nevertheless, using this wavefunction as a trial state, the variational energy is 
then given by 

h^V^ 1 Ze^' 
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E 



2m 
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Now, according to the variational principle, we must minimize the energy 
functional by varying with respect to the complex wavefunction, tpi, 

subject to the normalization condition, = 1. The latter can be imposed 

using a set of Lagrange multipliers, e^, i.e. 



E 



1 



0. 
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Following the variation,* one obtains the Hartree equations. 



2v72 



2m 



1 Ze' 



■4>i 



47re, 



-E 



/m2 



d'r\^P,{v')\ 



eiVi(r) . (9.4) 



Then according to the variational principle, amongst all possible trial functions 
tpi^ the set that minimizes the energy are determined by the effective potential. 




Douglas Rayner Hartree FRS 
1897-1958 

An English 
mathematician 
and physicist 
most famous for 
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analysis and its 
application to 
atomic physics. 
He entered 
St John's College Cambridge in 
1915 but World War I interrupted 
his studies and he joined a team 
studying anti-aircraft gunnery. He 
returned to Cambridge after the war 
and graduated in 1921 but, perhaps 
because of his interrupted studies, 
he only obtained a second class 
degree in Natural Sciences. In 1921, 
a visit by Niels Bohr to Cambridge 
inspired him to apply his knowledge 
of numerical analysis to the solution 
of differential equations for the 
calculation of atomic wavefunctions. 



^^W = i^E/'^V|^,(rOP-f^. 

Equation (9.4) has a simple interpretation: The first two terms relate to the 
nuclear potential experienced by the individual electrons, while the third term 
represents the electrostatic potential due to the other electrons. However, to 
simplify the procedure, it is useful to engineer the radial symmetry of the 
potential by replacing Ui{v) by its spherical average. 



Ui{v) ^ Ui{r) = 




®Note that, in applying the variation, the wavefunction ^* can be considered independent 
of - you might Uke to think why. 
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Finally, to relate the Lagrange multipliers, (which have the appearance 
of one-electron energies), to the total energy, we can multiply Eq. (9.4) by 



?/'j*(r) and integrate. 



2m 47reo r 



,2 



If we compare this expression with the variational state energy, we find that 
^-E^'-4^E / dVd^r\Ur')\^^mr)\'. (9-5) 

i i<j 

To summarize, if we wish to implement the central field approximation to 
determine the states of a multi-electron atom, we must follow the algorithm: 

1. Firstly, one makes an initial "guess" for a (common) central potential, 
U{r). As r ^ 0, screening becomes increasingly ineffective and we ex- 
pect U{r) ^0. As r — > oo, we anticipate that U{r) — > j^^^^r^j 
corresponding to perfect screening. So, as a starting point, we make 
take some smooth function U (r) interpolating between these limits. For 
this trial potential, we can solve (numerically) for the eigenstates of the 
single-particle Hamiltonian. We can then use these states as a plat- 
form to build the product wavefunction and in turn determine the self- 
consistent potentials, Ui{r). 

2. With these potentials, Ui{r), we can determine a new set of eigenstates 
for the set of Schrodinger equations. 



2m 47reo r 



3. An estimate for the ground state energy of an atom can be found by 
filling up the energy levels, starting from the lowest, and taking account 
of the exclusion principle. 

4. Using these wave functions, one can then make an improved estimate of 
the potentials Ui{ri) and return to step 2 iterating until convergence. 

Since the practical implemention of such an algorithm demands a large degree 
of computational flair, if you remain curious, you may find it useful to refer to 
the Mathematica code prepared by Ref. [4] where both the Hartree and the 
Hartree-Fock procedures (described below) are illustrated. 

> Info. An improvement to this procedure, known the Hartree-Fock method, 

takes account of exchange interactions. In order to do this, it is necessary to ensure 
that the wavefunction, including spin, is antisymmetric under interchange of any pair 
of electrons. This is achieved by introducing the Slater determinant. Writing the 
individual electron wavefunction for the i"^ electron as tjjkivi), where i = 1,2 - ■■ N and 
k is shorthand for the set of quantum numbers {nimems), the overall wavefunction 
is given by 



'/'i(ri) V'i(r2) V'i(r3) 

V'2(ri) V'2(r2) V'2(r3) 
V'3(ri) ?/'3(r2) tpsirz) 
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Note that each of the A''! terms in $ is a product of wavefunctions for each individual 
electron. The 1/\/~N'\ factor ensures the wavefunction is normalized. A determinant 
changes sign if any two columns are exchanged, corresponding to (say); this 

ensures that the wavefunction is antisymmetric under exchange of electrons i and j. 
Likewise, a determinant is zero if any two rows are identical; hence all the V'fcS must 
be different and the Pauli exclusion principle is satisfied.^ In this approximation, a 
variational analysis leads to the Hartree-Fock equations (exercise), 



£jl/'i(l') = 



2rn 



iirer. 



+ 



1 



47reo |r — r 



■ V'j(r)V'i(r')^n 



The first term in the last set of brackets translates to the ordinary Hartree contribution 
above and describes the influence of the charge density of the other electrons, while 
the second term describes the non-local exchange contribution, a manifestation of 
particle statistics. 

The outcome of such calculations is that the eigenfunctions cirG ^ cis for 
hydrogen, characterized by quantum numbers n, m^, with i < n, but that the 
states with different i for a given n are not degenerate, with the lower values 
of £ lying lower. This is because, for the higher £ values, the electrons tend 
to lie further from the nucleus on average, and are therefore more effectively 
screened. The states corresponding to a particular value of n are generally 
referred to as a shell, and those belonging to a particular pair of values of n, £ 
are usually referred to as a subshell. The energy levels are ordered as below 
(with the lowest lying on the left): 



Subshell name 


Is 


2s 


2p 


3s 


3p 


4s 


3d 


4p 


5s 


4d 


n = 


1 


2 


2 


3 


3 


4 


3 


4 


5 


4 


£ = 








1 





1 





2 


1 





2 


Degeneracy 


2 


2 


6 


2 


6 


2 


10 


6 


2 


10 


Cumulative 


2 


4 


10 


12 


18 


20 


30 


36 


38 


48 



Note that the values of Z corresponding to the noble gases, 2, 10, 18, 36, at 
which the ionization energy is unusually high, now emerge naturally from this 
filling order, corresponding to the numbers of electrons just before a new shell 
(n) is entered. There is a handy mnemonic to remember this filling order. By 
writing the subshells down as shown right, the order of states can be read off 
along diagonals from lower right to upper left, starting at the bottom. 

We can use this sequence of energy levels to predict the ground state 
electron configuration of atoms. We simply fill up the levels starting from 
the lowest, accounting for the exclusion principle, until the electrons are all 
accommodated (the aufbau principle). Here few examples: 



z 


Element 


Configuration 




loniz. Pot. (eV) 


i 


H 


(Is) 




13.6 


2 


He 


{Isf 




24.6 


3 


Li 


He {2s) 
He (2s)2 

He {2sf{2p) 
He (2s)2(2p)2 
He (2s)2(2p)3 
He (2s)2 (239)4 
He (2s)2(2p)'' 
He {2s)2(2p)6 




5.4 


4 


Be 




9.3 


5 


B 


Pl/2 


8.3 


6 


C 




11.3 


7 


N 




14.5 


8 





'P2 


13.6 


9 


F 


\P,/2 


17.4 


iO 


Ne 


So 


21.6 


11 


Na 


Ne (3s) 


'S,/2 


5.1 



tion 



^Note that for N = 
^[V'i(ri)V'2(r2) 



2, the determinant reduces to the familiar antisymmetric wavefunc- 

- V'2(l"l)V'l(l"2)]. 
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PERIODIC TABLE OF THE ELEMENTS 
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Li 



Na 



Rni.ilTlvn ATWIlt MASS (I 
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13 



Be 



Mg 



5 10.511 

-B 



|_i Metel Q Semimetal Qj Nonmeial 




lU Alkali moi4l 


^ CnaFi»j4l^» CiMKhl 










[j TfBnaikjr melali 


^ Noble 9 as 




[_j LanlhanidB 
□ Aclinirto 


SIMJDARD STATE {25 X. Id t'P-P} 



13 IIIA 14 IVA 15 VA 16 VIA 17 VIIA 



B 



:]I3 4 :V.ji 5 Vfi 6 VIB 7 VIIB 8 



Al 



Si 



O 



F 



CI 



He 

HELIUM 



Ne 



Ar 



19 39.0S 

K 



20 4CI.&S-8 

Ca 



Sc 



22 47My 

Ti 



23 M.342 

V 



24 SLSae 

Cr 



25 W.^M 

Mil 



26 S5.S45 

Fe 



27 SS.Si3 28 3S,693 

Co Ni 



29 ^.fAS 

Cu 



30 

Zn 



31 

Ga 



32 ra.fr 

Ge 



33 J4.9S2 

As 



34 ?6,9i 

Se 



35 7S.S« 

Br 



36 &3.Si> 

Kr 



37 85.4ee 

Rb 



Sr 



39 aB.90B 

Y 



Zr 

zinccNim 



^41 9ij.SD£ 

Nb 

NIOBIUM 



41 flS94 

Mo 
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Ru 

nUTHENilUU 



*5 1&2.91 46 IM-dS 

Rh Pd 



Ag 



Cd 

-CAOMiOM 



In 

IMO-:l>^ 



50 118.T 

Sn 



51 i!i7e 

Sb 



52 12769 

Te 



I 

IODINE 



54 131.Sfl 

Xe 



55 isa 91 

Cs 



Ba 



57-71 
l.a 

.iiiitli^nEdi 



72 1»8,49 

Hf 



73 ISO, 95 

Ta 



74 183,8^ 

w 



Re 



76 130,2; 

Os 

OSMIUW 



77 192 £2 Its i95,M 

Ir Pt 



Au 



Hg 



SI 204,3? 

Tl 



Pb 



83 M,9 

Bi 



S4 (£09> 

Po 

POLONHJMi 



At 



86 1222) 

Rn 



Fr 
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Ra 
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1(19 (26aj 110 (3B1) 
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Ho 


6S 167,26 

Er 

iRBIIJtH 


69 i&e,93 

Tm 


70 173.oa 

Yb 

■iTiERSigAi 


71 174.97 
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Figure 9.3: Periodic table of elements. 



Since it is generally the outermost electrons which are of most interest, con- 
tributing to chemical activity or optical spectra, one often omits the inner 
closed shells, and just writes O as (2p)^, for example. However, the configura- 
tion is not always correctly predicted, especially in the heavier elements, where 
levels may be close together. It may be favourable to promote one or even 
two electrons one level above that expected in this simple picture, in order to 
achieve a filled shell. For example, Cu {Z = 29) would be expected to have con- 
figuration ■ ■ ■ (4s)^(3d)^, and actually has configuration ■ ■ ■ (4s)-'^(3d)^°. There 
are several similar examples in the transition elements where the d subshells 
are being filled, and many among the lanthanides (rare earths) and actinides 
where / subshells are being filled. 

> Info. Since the assignment of an electron configuration requires only the 
enumeration of the values of n and £ for all electrons, but not those of me and rus, 
each configuration will be accompanied by a degeneracy g. If Vnt denotes the number 
of electrons occupying a given level E^/, and 6e = 2 x {2i + 1) is the degeneracy of 
that level, there are 

dni = ^ (9.6) 

ways of distributing the Vni, electrons among the 5t individual states. The total 
degeneracy, is then obtained from the product. 

This scheme provides a basis to understand the periodic table of el- 
ements (see Fig. 9.3). We would expect that elements which have similar 
electron configurations in their outermost shells (such as Li, Na, K, Rb, Cs, Fr 
which all have (ns)^ or F, CI, Br, I, which all have [np)^) would have similar 
chemical properties, such as valency, since it is the unpaired outer electrons 
which especially participate in chemical bonding. Therefore, if one arranges 
the atoms in order of increasing atomic number Z (which equals the number 
of electrons in the atom), periodic behaviour is seen whenever a new subshell 
of a given £ is filled. 
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9.3 Coupling schemes 

The procedure outlined above allows us to predict the occupation of subshells 
in an atomic ground state. This is not in general sufficient to specify the ground 
state fully. If there are several electrons in a partially filled subshell, then their 
spins and orbital angular momenta can combine in several different ways, to 
give different values of total angular momentum, with different energies. In 
order to deal with this problem, it is necessary to consider the spin-orbit 
interaction as well as the residual Coulomb interaction between the outer 
electrons. 

Schematically we can write the Hamiltonian for this system as follows: 

i<j ^■^ i i 

V ^ ' ^ V ' 

Hi H2 

where Hq includes the kinetic energy and central field terms, Hi is the residual 

Coulomb interaction, and (with £,i{ri) = 2^^^2^ I: {drV (r))) H2 is the spin-orbit 
interaction. We can then consider two possible scenarios: 

Hi ^ H2: This tends to apply in the case of light atoms. In this situation, 
one considers first the eigenstates of Hq + Hi, and then treats H2 as a 
perturbation. This leads to a scheme called LS (or Russell-Saunders) 
coupling. 

H2 ^ Hi'. This can apply in very heavy atoms, or in heavily ionized light 
atoms, in which the electrons are moving at higher velocities and rela- 
tivistic effects such as the spin-orbit interaction are more important. In 
this scheme called jj coupling applies. 

It is important to emphasise that both of these scenarios represent approxima- 
tions; real atoms do not always conform to the comparatively simple picture 
which emerges from these schemes, which we now discuss in detail. 

9.3.1 LS coupling scheme 

In this approximation, we start by considering the eigenstates of Hq + Hi. We 
note that this Hamiltonian must commute with the total angular momentum 
(because of invariance under rotations in space) , and also clearly commutes 
with the total spin S^. It also commutes with the total orbital angular momen- 
tum L^, since Hi only involves internal interactions, and must therefore be 
invariant under global rotation of all the electrons. Therefore the energy levels 
can be characterised by the corresponding total angular momentum quantum 
numbers L, S, J. Their ordering in energy is given by Hund's rules: 

1. Combine the spins of the electrons to obtain possible values of total spin 
S. The largest permitted value of S lies lowest in energy. 

2. For each value of 5, find the possible values of total angular momentum 
L. The largest value of L lies lowest in energy. 

3. Couple the values of L and S to obtain the values of J (hence the name 
of the scheme). If the subshell is less than half full, the smallest value 
of J lies lowest; otherwise, the largest value of J lies lowest. 
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In deciding on the permitted values of L and S, in addition to applying the 
usual rules for adding angular momenta, one also has to ensure that the exclu- 
sion principle is respected, as we will see later when considering some examples. 

These rules are empirical; there are exceptions, especially to the L and J 
rules (2 and 3). Nevertheless, Hund's rules are a useful guide, and we should 
try to understand their physical origin. 

1. Maximising S makes the spin wavefunction as symmetric as possible. 

This tends to make the spatial wavefunction antisymmetric, and hence 
reduces the Coulomb repulsion, as we saw when discussing the exchange 
interactions in Helium. 

2. Maximising L also tends to keep the electrons apart. This is less obvious, 
though a simple classical picture of electrons rotating round the nucleus 
in the same or different senses makes it at least plausible. 

3. The separation of energies for states of different J arises from treating 

the spin-orbit term H2 as a perturbation (fine structure). It can be 
shown (using the Wigner-Eckart theorem - beyond the scope of these 
lectures) that 

{\J,mj,L,S\^ii{ri)Li ■ Si\J,mj,L,S) 

i 

= C{L, S){J, mj, L, S\t . S| J, mj, L, S) 

= ^-^^W + 1) - L{L + 1) - S{S + 1)] , (9.7) 

where the matrix element C(L, S) depends on the total L and S values. 
Since one may show that the sign of C,{L, S) changes according to the 
whether the subshell is more or less than half-filled, the third Hund's 
rule is established. 

To understand the application of LS coupling, it is best to work through 
some examples. Starting with the simplest multi-electron atom, helium, the 
ground state has an electron configuration (Is)^, and must therefore have 
L — S — J — {). In fact, for any completely filled subshell, we have L — S — Q 
and hence J = 0, since the total and mg must equal zero if all substates 
are occupied. Consider now an excited state of helium, e.g. (ls)^(2p)^, in 
which one electron has been excited to the 2p level. We can now have S = \ 
or S" = 0, with the S = 1 state lying lower in energy according to Hund's 
rules. Combining the orbital angular momenta of the electrons yields L — 1 
and thus, with S — J — 1, while with S — 1, J = 0, 1, 2 with J — Q lying 
lowest in energy. 

Once again, as with the hydrogen-like states, we may index the states of 
multi-electron atoms by spectroscopic term notation, ^'^+^Lj. The superscript 
2S + \ gives the multiplicity of J values into which the level is split by the 
spin-orbit interaction; the L value is represented by a capital letter, 5, P, D, 
etc., and J is represented by its numerical value. Thus, for the (ls)^(2p)^ 
state of helium, there are four possible states, with terms: 

=^Po 'Pi 'P2 'Pi , 

where the three ^P states are separated by the spin-orbit interaction, and the 
singlet ^P state lies much higher in energy owing to the Coulomb interaction. 
The separations between the 'P2 and 'Pi and the 'Pi and 'Pq should be in the 
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ratio 2:1. This is an example of the Lande interval rule, which states that 
the separation between a pair of adjacent levels in a fine structure multiplet is 
proportional to the larger of the two J values involved. This is easily shown 
using Eq. (9.7) - the separation in energy between states J and J — 1 is 

oc J(J+ 1) - (J- 1)J= 2J. 

Actually in the case of helium the situation is a bit more complicated, because 
it turns out that the spin-orbit interaction between different electrons makes 
a non-negligible additional contribution to the fine structure. Other excited 
states of helium, of the form (Is)^(n^)^, can be handled similarly, and again 
separate into singlet and triplet states. 

O Exercise. For the case of boron, with the electron configuration (ls)^(2s)^(2p), 
use Hund's rules to show that the ground state is 

We next consider the case of carbon, which has ground state electron 
configuration (ls)^(2s)^(2p)^. This introduces a further complication; we now 
have two identical electrons in the same unfilled subshell, and we need to 
ensure that their wavefunction is antisymmetric with respect to electron ex- 
change. The total spin can either be the singlet 5 = state, which has an 
antisymmetric wavefunction -^[\ Ti) | i2) — | ii) <S) | T2)] , or one of the triplet 
S =1 states, which are symmetric, ^[| ti) ® | 12) + Ui) ® | T2)], I Ti) ® I T2) 
or I ii) <8) I 12)- We must therefore choose values of L with the appropriate 
symmetry to partner each value of S. To form an antisymmetric state, the two 
electrons must have different values of m^, so the possibilities are as shown 
right. Inspecting the values of tjil we can deduce that L = 1.^" By contrast, 
to form a symmetric total angular momentum state, the two electrons may 
have any values of m^, leading to the possibilities shown right. Inspecting the 
values of mi we infer that L = 2 or 0. 

We must therefore take S — 1 with L — 1 and S — with L = 2 or 0. 
Finally, to account for the fine structure, we note that the states with S — 1 
and L — 1 can be combined into a single J = state, three J — 1 states, and 
five J = 2 states leading to the terms ^Pq, ^Pi, and •^P2 respectively. Similarly 
the S — 0, L = 2 state can be combined to give five J = 2 states, ^D2, while the 
S — L — Q state gives the single J = state, ^Sq. Altogther we recover the 
1 + 3 + 5 + 5 + 1 = 15 possible states (cf. Eq. (9.6) with the ordering in energy 
given by Hund's rules (shown to the right). The experimental energy values 
are given using the conventional spectroscopic units of inverse wavelength. 
Note that the Lande interval rule is approximately obeyed by the fine structure 
triplet, and that the separation between L and S values caused by the electron- 
electron repulsion is much greater than the spin-orbit effect. 

In an excited state of carbon, e.g. (2p)^(3p)^, the electrons are no longer 
equivalent, because they have different radial wavef unctions. So now one can 
combine any of S* = 0, 1 with any of L = 0, 1, 2, yielding the following terms 
(in order of increasing energy, according to Hund's rules): 

■^01,2,3 ^Po,i,2 ^Si ^Pi ^So . 

For nitrogen, the electron configuration is given by (ls)^(2s)^(2p)3. The 
maximal value of spin is 5 = 3/2 while L can take values 3, 2, 1 and 0. Since 

'^"This result would also be apparent if we recall the that angular momentum states are 
eigenstates of the parity operator with eigenvalue (—1)^. Since there axe just two electrons, 
this result shows that both the L = and L = 2 wavefunction must be symmetric under 
exchange. 
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the spin wavefunction (being maximal) is symmetric, the spatial wavefunction 
must be completely antisymmetric. This demands that all three states with 
= 1,0,-1 must be involved. We must therefore have L = 0, leading to 
J = 3/2 and the term, 

> Exercise. Construct the L = state involving the addition of three £ = 1 
angular momentum states. Hint: make use of the total antisymmetry condition. 



's % 



(Pill''' 



•->2 — -1 



As a final example, let us consider the ground state of oxygen, which has 
electron configuration {2p)'^. Although there are four electrons in the (2p) 
subshell, the maximum value of 5 = 1. This is because there are only three 
available values of = ±1,0, and therefore one of these must contain two 
electrons with opposite spins. Therefore, the maximum value of mg = 1, 
achieved by having electrons with = +^ in both the other states. By 
pursuing this argument, it is quite easy to see that the allowed values of L, S 
and J are the same as for carbon (2p)^. This is in fact a general result - the 
allowed quantum numbers for a subshell with n electrons are the same as for 
that of a subshell with n "holes" . Therefore, the energy levels for the oxygen 
ground state configuration are the same as for carbon, except that the fine 
structure multiplet is inverted, in accordance with Hund's third rule. 

9.3.2 jj coupling scheme 

When relativitic effects take precedence over electron interaction effects, we 
must start by considering the eigenstates of Hq + H2 = Hq + ^ • ^i(ri)Li ■ S^. 
These must be eigenstates of as before, because of the overall rotational 
invariance, and also of J? for each electron. Therefore, in this case, the cou- 
pling procedure is to find the allowed j values of individual electrons, whose 
energies will be separated by the spin-orbit interaction. Then these individual 
j values are combined to find the allowed values of total J. The effect of the 
residual Coulomb interaction will be to split the J values for a given set of js. 
Sadly, in this case, there are no simple rules to parallel those of Hund. 

As an example, consider a configuration (np)^ in the jj coupling scheme, to 
be compared with the example of carbon which we studied in the LS scheme. 
Combining s = 1/2 with i = 1, each electron can have j = 1/2 or 3/2. If 
the electrons have the same j value, they are equivalent, so we have to take 
care of the symmetry of the wavefunction. We therefore have the following 
possibilities: 

^ Ji = j2 = 3/2 => J = 3, 2, 1,0, of which J = 2, are antisymmetric. 

> Ji = j2 = 1/2 => J = 1, 0, of which J = is antisymmetric. 

> Ji = l/2, i2 = 3/2^ J = 2,l. 

In jj coupling, the term is written {ji,j2).j, so we have the following terms in 
our example: 



(1/2,1/2)0 (3/2,1/2)1(3/2,1/2)2 



(3/2,3/2)2 (3/2,3/2)0 



in order of increasing energy. Note that both LS and jj coupling give the same 
values of J (in this case, two states with J = 0, two with J = 2 and one with 
J = 1) and in the same order. However, the pattern of levels is different; in 
LS coupling we found a triplet and two singlets, while in this ideal jj scenario, 
we have two doublets and a singlet. The sets of states in the two coupling 



Level scheme of the carbon 
atom (ls)2(2s)2(2p)2. Draw- 
ing is not to scale. On the 
left the energy is shown with- 
out any two-particle interaction. 
The electron-electron interaction 
leads to a three-fold energy split- 
ting with L and S remaining 
good quantum numbers. Spin- 
orbit coupling leads to a further 
splitting of the states with J re- 
maining a good quantum num- 
ber. Finally on the right, the lev- 
els show Zeeman splittings in an 
external magnetic field. In this 
case, the full set of 15 levels be- 
come non-degenerate. 
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schemes must be expressible as linear combinations of one another, and the 
physical states for a real atom are likely to differ from either approximation. 

In fact, this idealized form of jj coupling is not seen in the heaviest such 
atom in the periodic table, lead (6p)^. However, it is seen in some highly 
ionized states, for example in Cr^^+, which has the same electron configu- 
ration (2p)^ as carbon, but where, because of the larger unscreened charge 
on the nucleus, the electrons are moving more relativistically, enhancing the 
spin-orbit effect. However, a classic example of the transition from LS to jj 
coupling is seen in the series C-Si-Ge-Sn-Pb in the excited states (2p)(3s), 
(3p)(4s), ■ ■ ■ (6p)(7s) (see figure right). Here, the electrons are not in the same 
subshell, so their wavefunctions overlap less, and the Coulomb repulsion is 
reduced compared to the spin-orbit interaction. Analysing this situation in 
the LS coupling approximation, one expects a triplet and a singlet: 

'P0,l,2 'Pi , 

while in the jj scheme one expects two doublets: 



Pure LS 



Pure JJ 



(1/2,1/2)0- 



:i/2, 3/2)2, 



1 • 



Experimentally, C and Si conform to the LS expectation and Pb to the jj 
scheme, while Ge and Sn are intermediate. 




Si Ge Sn 



9.4 Atomic spectra 

Atomic spectra result from transitions between different electronic states of 
an atom via emission or absorption of photons. In emission spectra, an 
atom is excited by some means (e.g. thermally through collisions, or by an 
electric discharge) , and one observes discrete spectral lines in the light emitted 
as the atoms relax. In absorption spectra, one illuminates atoms using a 
broad waveband source, and observes dark absorption lines in the spectrum 
of transmitted light. Of course the atoms excited in this process subsequently 
decay by emitting photons in random directions; by observing in directions 
away from the incident light this fluorescence radiation may be studied. The 
observation of these spectral lines is an important way of probing the atomic 
energy levels experimentally. In the case of optical spectra and the nearby 
wavebands, the excited states responsible generally involve the excitation of 
a single electron from the ground state to some higher level. In some cases, 
it may simply involve a different coupling of the angular momenta within the 
same electron configuration. These are the kinds of excitations which we are 
about to discuss. However, other types of excitations are also possible. For 
example. X-ray emission occurs when an electron has been removed from one 
of the innermost shells of a heavy atom; as electrons cascade down to fill the 
hole, high energy photons may be emitted. 

The basic theory governing stimulated emission and absorption, and spon- 
taneous emission of photons will be outlined in detail when we study radiative 
transitions in chapter 13. Here we must anticipate some of the basic con- 
clusions of that study. In the electric dipole approximation, the rate of 
transitions is proportional to the square of the matrix element of the electric 
dipole operator between the initial and final states, |(?/^f|d|?/'i)p. In addition, 
the rate of spontaneous transitions is proportional to cu^, where a; = \E{ — Ei\ 
denotes the energy separation between the states. 

The form of the dipole operator, d means that the matrix elements may 
vanish identically. This leads to a set of selection rules defining which tran- 
sitions are allowed. Here we consider the simplest case of a single electron, but 
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the principles can be generalized. Referring to chapter 13 for a more detailed 
discussion, one finds that, for a transition to take place: 

> Parity must change; 

> AJ = ±1, (but ^ is not allowed) and AMj = ±1, . 

Atomic states are always eigenstates of parity and of total angular momentum, 
J, so these selection rules can be regarded as absolutely valid in electric dipole 
transitions. It should be emphasized again, though, that the electric dipole 
approximation is an approximation, and higher order processes may occur, 
albeit at a slower rate, and have their own selection rules. 

In specific coupling schemes, further selection rules may apply. In the case 
of ideal LS coupling, we also require: 

> AS = and AMs = 0; 

> AL = ±1, (but ^ is not allowed) and AMl = ±1, 0; 

> and A£i = ±1 if only electron i is involved in the transition. 

In LS coupling, the states are eigenstates of total spin; since the dipole operator 
does not operate on the spin part of the wavefunction, the rules on AS and 
AMs follow straightforwardly. This, and the absolute rules relating to J, 
imply the rules for L and Ml. The rule for A£i follows from the parity change 
rule, since the parity of the atom is the product of the parities of the separate 
electron wavef unctions, given by ( — 1)^'. However, since LS coupling is only 
an approximation, these rules should themselves be regarded as approximate. 
With this preparation, we now turn to the consequences of the selection rules 
on the atomic spectra of atoms. 



9.4.1 Single electron atoms 

In this context, "single electron atoms" refer to atoms whose ground state 
consists of a single electron in an s level, outside closed shells; it is this electron 
which is active in optical spectroscopy. Our discussion therefore encompasses 
the alkali metals, such as sodium, and also hydrogen. We take sodium, 

whose ground state configuration is (3s) ^, as our example: 

3 04 Spin-orbit QQ2'\ eV ^ 

> The ground state has term ^Si/2- The excited states are all doublets splitting 
with J = L± 1/2, except for the s states, which are obviously restricted 
to J= 1/2. 

> The parity is given by ( — 1)^', so the allowed transitions involve changes ^ 
in £ by ±1 unit, i.e. s ^ p, p ^ d, d ^ f, etc. Changes of more than S 

one unit in £ would fall foul of the A J rule. i <S ,i ko.597 nm 

t> The s ^ p transitions are all doublets. All the doublets starting or 
ending on a given p state have the same spacing in energy. The transition 
3s <-> 3p gives rise to the familiar yellow sodium "D-lines" at 589 nm .^.i^ 
(see right). 

[> The p ^ d transitions involve two doublets, ^Pi/2,3/2 ^03/2,5/2- 
However, the ^^05/2 transition is forbidden by the AJ rule, so 

the line is actually a triplet. In practice, the spin-orbit interaction falls 
quite rapidly with increasing £ (and with increasing n) as the effect of 
screening increases, so that the effect of the ^03/2,5/2 splitting may not 
be resolved experimentally. 
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> As n increases, the energy levels approach (from below) those for hy- 
drogen, because the nuclear charge is increasingly effectively screened 
by the inner electrons. This happens sooner for the higher £ values, for 
which the electron tends to lie further out from the nucleus. 

i> In an absorption spectrum, the atoms will start from the ground state, 
so only the 3s — > np lines will be seen. In emission, the atoms are excited 
into essentially all their excited levels, so many more lines will be seen 
in the spectrum. 

The comments above for sodium also apply for hydrogen, except that, in this 
case, (2s, 2p), (3s, 3p, 3d), etc. are degenerate. One consequence is that the 2s 
state in hydrogen is metastable - it cannot decay to the only lower lying level 
(Is) by an electric dipole transition. In fact its favoured spontaneous decay 
is by emission of two photons; a process which is described by second-order 
perturbation theory. In practice, hydrogen atoms in a 2s state are more likely 
to deexcite through collision processes. During an atomic collision, the atoms 
are subject to strong electric fields, and we know from our discussion of the 
Stark effect that this will mix the 2s and 2p states, and decay from 2p to Is 
is readily possible. 

9.4.2 Helium and alkali earths 

We next discuss atoms whose ground state consists of two electrons in an s 
level. Our discussion therefore covers helium (Is)^, and the alkali earths: 
beryllium (2s) ^, magnesium (3s) ^, calcium (4s) ^, etc. We start with helium. 

The ground state has term ^Sg. The excited states are of the form 
(ls)(n£) (the energy required to excite both of the Is electrons to higher 
states is greater than the first ionization energy, and therefore these 
form discrete states within a continuum of ionized He^+e" states). The 
excited states can have S' = or 5 = 1, with S — 1 lying lower in energy 
(Hund). 

> The LS coupling approximation is pretty good for helium, so the AS" = 
selection rule implies that the S = and S = 1 states form completely 
independent systems as far as spectroscopy is concerned. 

t> The lines in the S = system are all singlets. They can be observed 
in emission, and those starting from the ground state can be seen in 
absorption. 

> The lines in the S — 1 system are all multiplets. They can be observed 
in emission only. Transitions of the form ^Si ^^P2.i,o are observed as 
triplets, spaced according to the Lande interval rule. Transitions of the 
form ^P2,i,o •*-^^D3,2,i are observed as sextuplets, as is easily seen by 
application of the AJ — ±1, rule. Actually, as mentioned above, the 
fine structure is a little more subtle in the case of helium. 

The alkali earths follow the same principles. In the case of calcium, the 
triplet 4p state is the lowest lying triplet state, and therefore metastable. In 
fact a faint emission line corresponding to the ^Pi — >^So decay to the ground 
state may be observed; this violates the AS — rule, and indicates that the LS 
coupling approximation is not so good in this case. A more extreme example is 
seen in Mercury, ground state (6s)^(5d)^°. Excited states involving promotion 
of one of the 6s electrons to a higher level can be treated just like the alkali 
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earths. In this case the "forbidden" ^Pi — >^ 5o is actually a prominent feature 
of the emission spectrum in the visible, implying a significant breakdown of 
the LS approximation. 



9.4.3 Multi-electron atoms 

Similar principles can be used to make sense of the spectra of more complicated 
atoms, though unsurprisingly their structure is more complex. For example, 
carbon, with ground state (2s)^(2p)^, corresponds to terms ^Po,i,2, and 
^So as discussed above. The excited states are of the form (2s)^(2j3)^(ni!)^, 
and can be separated into singlets and triplets, and in addition excitations 
of the form (2s)^(2p)^ can arise. Nitrogen, with ground state (2s)^(2p)^, has 
three unpaired electrons, so the ground state and excited states form doublets 
(S = 1/2) and quartets (5 = 3/2) with correspondingly complex fine structure 
to the spectral lines. 



9.5 Zeeman effect 
9.5.1 Single-electron atoms 

Before leaving this chapter on atomic structure, we are now in a position to 
revisit the question of how atomic spectra are influenced by a magnetic field. 
To orient our discussion, let us begin with the study of hydrogen- like atoms 
involving just a single electron. In a magnetic field, the Hamiltonian of such 
a system is described by ^ = + -^rei. + -ffzeemam where denotes the 
non-relativistic Hamiltonian for the atom, Hrei incorporates the relativistic 
corrections considered earlier in the chapter, and 

eB - 

-ffzeeman — ~n [Lz + 2S'z) , 

2mc 

denotes the Zeeman energy associated with the coupling of the spin and orbital 
angular momentum degrees of freedom to the magnetic field. Here, since 
we are dealing with confined electrons, we have neglected the diamagnetic 
contribution to the Hamiltonian. Depending on the scale of the magnetic field, 
the spin-orbit term in Hrei. or the Zeeman term may dominate the spectrum 
of the atom. 

Previously, we have seen that, to leading order, the relativistic corrections 
lead to a fine-structure energy shift of 



AET- = -mc' — 



"■■^ 2 \ n J \4 3 + 1/2^ 

for states |n,j — l± l/2,mj,i). For weak magnetic fields, we can also treat 
the Zeeman energy in the framework of perturbation theory. Here, although 
states with common j values (such as 25i/2 a^i^d 2Pxj2) are degenerate, the 
two spatial wavefunctions have different parity (£ = and 1 in this case), 
and the off-diagonal matrix element of .ffzeeman coupling these states vanishes. 
We may therefore avoid using degenerate perturbation theory. Making use of 
the relation (exercise - refer back to the discussion of the addition of angular 
momenta and spin in section 6.4.2), 

{n,j = i ± l/2,mj,£\Sz\n,j = i±l/2,mj,i) = , 

we obtain the following expression for the first order energy shift, 

AEfZT^r;^.,, = i ± 1/2, mj,i) = MBSm,- (l ± ^) , 
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Figure 9.4: The well known doublet which is responsible for the bright yellow light 
from a sodium lamp may be used to demonstrate several of the influences which 
cause splitting of the emission lines of atomic spectra. The transition which gives 
rise to the doublet is from the 3p to the 3s level. The fact that the 3s state is lower 
than the 3p state is a good example of the dependence of atomic energy levels on 
orbital angular momentum. The 3s electron penetrates the Is shell more and is less 
effectively shielded than the 3p electron, so the 3s level is lower. The fact that there 
is a doublet shows the smaller dependence of the atomic energy levels on the total 
angular momentum. The 3p level is split into states with total angular momentum 
J = 3/2 and J = 1/2 by the spin-orbit interaction. In the presence of an external 
magnetic field, these levels are further split by the magnetic dipole energy, showing 
dependence of the energies on the ^-component of the total angular momentum. 



where /^b denotes the Bohr magneton. Therefore, we see that all degenerate 
levels are split due to the magnetic field. In contrast to the "normal" Zeeman 
effect, the magnitude of the splitting depends on L 

t> Info. If the field is strong, the Zeeman energy becomes large in comparison 
with the spin-orbit contribution. In this case, we must work with the basis states 
\n^£,me,ms) = \n,£,me) ® Im^) in which both Hq and ifzeeman are diagonal. Within 
first order of perturbation theory, one then finds that (exercise) 

in/ Za\'^ /3 n nrufnig \ 

the first term arising from the Zeeman energy and the remaining terms from ifrei.- At 
intermediate values of the field, we have to apply degenerate perturbation theory to 
the states involving the linear combination of | n, j = £± 1/2, mj,£). Such a calculation 
reaches beyond the scope of these lectures and, for details, we refer to the literature 
(see, e.g., Ref. [6]. Let us instead consider what happens in multi-electron atoms. 



9.5.2 Multi-electron atoms 

For a multi-electron atom in a weak magnetic field, the appropriate unper- 
turbed states are given by \J,Mj,L,S), where J, L, S refer to the total 
angular momenta. To determine the Zeeman energy shift, we need to deter- 
mine the matrix element of Sz- To do so, we can make use of the following 
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argument. Since 2L ■ S = — — S^, this operator is diagonal in the basis of 
states, \ J,Mj,L,S). Therefore, the matrix element of the operator (exercise, 
hint: recall that [Si,Sj] = ihcijkSk and [Li,Sk] = 0), 

-ihS X L = S(L ■ S) - (L ■ S)S 



must vanish. Moreover, from the identity [L 
matrix element of the vector product. 



S,J] = 0, it follows that the 



-ih{S X L) X J = S X J(L ■ S) - (L ■ S)S X J , 

must also vanish. If we expand the left hand side, we thus find that the matrix 
element of 

(S X L) X J = L(S ■ J) - S(L ■ J) ^=="^ J(S ■ J) - SJ2 , 

also vanishes. Therefore, it follows that (SJ^) = (J(S ■ J)), where the expec- 
tation value is taken over the basis states. Then, with S ■ J = ^(J^ + — L^), 
we have that 



{Jz 



J{J + I) + S{S + I) - L{L + I) 
2J(J + 1) 



As a result, we can deduce that, at first order in perturbation theory, the 
energy shift arising from the Zeeman term is given by 



where 



gj = i 



J{J + 1) + S{S+1)-L{L + 1) 
2J(J + 1) 



denotes the effective Lande g- factor, which lies between 1 and 2. Note that, 
in the special case of hydrogen, where 5" = 1/2 and J = L ± 1/2, we recover 
our previous result. The predicted Zeeman splitting for sodium is shown in 
figure 9.4. 



> Info. In the strong field limit, where the influence of Zeeman term domi- 
nates, the appropriate basis states are set by \L, M^, S, Ms), in which the operators 

Lz, S , Sz, and ifzeeman 

form 



L^, Lz, S'^, 5*2, and -ffzEEman are diagonal. In this case, the energy splitting takes the 



where the second term arises from the spin-orbit interaction. 



I Weak external J = |_ 4- 3 
B nnagnelic field , 




g Strang e>rternal 




itragnetic field 



Orbital and spin 
angular rrwmenta 
irtdep&ndantly 
precess about 
magnetic field 
direction J quanturri 
nunit>ei no longer 
a constant of the 
motion. 



Figure 9.5: In the weak field 
case, the vector model (top) im- 
plies that the coupling of the or- 
bital angular momentum L to 
the spin angular momentum S 
is stronger than their coupling 
to the external field. In this 
case where spin-orbit coupHng is 
dominant, they can be visual- 
ized as combining to form a to- 
tal angular momentum J which 
then precesses about the mag- 
netic field direction. In the 
strong field case, S and L cou- 
ple more strongly to the exter- 
nal magnetic field than to each 
other, and can be visualized as 
independently precessing about 
the external field direction. 
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